Anomalous transport and quantum-classical correspondence 
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We present evidence that anomalous transport in the classical standard map results in strong 
enhancement of fluctuations in the localization length of quasienergy states in the corresponding 
quantum dynamics. This generic effect occurs even far from the semiclassical limit and reflects 
the interplay of local and global quantum suppression mechanisms of classically chaotic dynamics. 
Possible experimental scenarios are also discussed. 
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It is generally accepted that quantum mechanics sup- 
presses chaotic classical motion. Numerous studies have 
identified mechanisms for suppression which, for our pur- 
poses, fall into two broad classes. One class is exemplified 
by 'dynamical localization' where quantum eigenstates 
are localized in an action variable like momentum despite 
the deterministic diffusion exhibited in the limiting classi- 
cal dynamics Q . This effect is analogous to the insulator 
phase or Anderson localization in tight-binding models 
H and, in one-dimension, the localized wavefunction has 
a characteristic exponential form. An associated scale is 
the localization length £ which is related to the classi- 
cal diffusion constant, an important fact for the work re- 
ported here. There also exists a demarcating timescale t* 
beyond which quantum and classical dynamics deviate. 
Despite a few counterexamples, dynamical localization 
provides a 'global' mechanism for quantum suppression. 

The second class can be motivated by the fast devi- 
ation of wavepacket dynamics from the classical motion 
even in the semiclassical limit. This effect is contained 
in the logarithmic break time t% = In (I /h)/T for the 
breakdown of quantum-classical correspondence Q| , rela- 
tive to a characteristic Lyapunov exponent T and action 
I. However, the quantum dynamics retains features of 
the classical dynamics for times beyond this estimate. 
The phenomenon of 'scarring' |^|,|| refers to quantum co- 
herences associated with local, unstable and marginally 
stable, classical invariant structures. Though there are 
still many open questions, this effect has been numeri- 
cally and experimentally observed in a wide variety of 
strongly coupled quantum systems. 

In this letter, we present evidence for another exam- 
ple of 'classical persistence' related to anomalous diffu- 
sion and accelerator modes in chaotic dynamics. Re- 
cently, there has been considerable activity in the area 
of anomalous classical transport |||7],|8| . What is of par- 
ticular relevance to the present work is that even in an 
established paradigm like the standard map, described 
by q n +i = q n + Pn+i and p n +i = p n + Ksmq n , some 
surprising new results were reported. 



In the limit of large K, it is well-known that the classi- 
cal standard map dynamics is diffusive in the action vari- 
able p with diffusion constant < p 2 > jt = K 2 /2 = D q i. 
What is also established is that away from this limit, the 
diffusion constant exhibits peaks with changing K and 
varies in time j|] . Recent results [^|J^] with improved pre- 
cision and detail show that there are values of K where 
the effective diffusion constant (over some time) can be 
different from D q i by many orders of magnitude. Fig- 
ure [l] shows these windows of anomalous diffusion which 
are not very narrow and exhibit sub-structure in each of 
the peaks. 
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FIG. 1. Diffusion coefficient D for the standard map, nor- 
malized to Dqi = K 2 /2, as a function of K. The deviations 
from the quasilinear prediction span a wide window of K and, 
as seen from the inset, can be large at select values of K, even 
far from the bifurcation point K — 2n. 

It was shown that there is no diffusion at all and that 
the real random walk process corresponds to a Levy-like 
wandering characterized by a transport exponent /i given 
by < p 2 >« where fj, > 1. /.i varies with K and has 
the 'normal' value fi — 1 only at special values of K. 
Superdiffusion occurs for fi > 1 leading to anomalous 
growth in momentum. It was explicitly shown Q that 
this could result from an island hierarchy with a peculiar 
topological structure near the 'accelerator' modes that 
appear at special values of K. Near these values, classical 
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trajectories stick to the boundaries of these islands which, 
in turn, leads to flights of arbitrary length. The net result 
is strong intermittency and superdiffusion. 

Another signature of anomalous transport results from 
the fact that a classical trajectory originating in a re- 
gion of phase space recurs infinitely often. From the 
set of recurrence times {r,}, j = 0, 1, • • •, Poincare cy- 
cles {tj} = {tj + i — Tj} and an associated probability 
distribution function R(t) can be constructed. In the 
case of 'perfect' mixing and normal diffusion, R(t) is 
strictly Poissonian whereas it exhibits power-like asymp- 
totics for anomalous kinetics. This classical characteris- 
tic of anomalous transport is illustrated in Fig. 0, and 
is important in the corresponding quantum dynamics as 
well. 




FIG. 2. Distribution of Poincare cycles. The inset clearly 
displays the power-law tail of the distribution. 

We raise several questions, arising from the classical 
anomalous diffusion, in the corresponding quantum dy- 
namics of the standard map or, equivalently, the 'delta- 
kicked rotor', H = p 2 /2 + K cos q ^2 n S(t — n) . The scale 
length £ for dynamical localization is related to the clas- 
sical diffusion constant and even follows the variations 
in D with K Jll]]. Thus, we anticipate that the strong 
enhancement in D resulting from anomalous diffusion 
should also be reflected in £. The nature of dynamical 
localization in these special windows is a related ques- 
tion. Further, as small local structures in the classical 
phase space lead to the superdiffusion, these issues di- 
rectly pertain to the interplay of local and global aspects 
of quantum suppression. This idea is reinforced by the 
enhanced diffusion that ought to be visible in the quan- 
tum dynamics even far from the semiclassical limit, where 
scarring is significant. Other authors have considered 
quantum dynamics in the presence of anomalous trans- 
port P,E0| though in a different regime. We concentrate 
exclusively on the situation where the quantization scale 
2ir?i is much larger than the size of the islands. Thus, 
any coherences associated with the islands would be due 
to 'scarred' quantum states. 

Consider the evolution of < p 2 > with time t (mea- 
sured as the number of kicks) in the quantum dynamics 
for two values of the classical stochasticity parameter K. 



The computation uses fast-Fourier transforms to evolve 
a plane wave initial state under repeated application of 
the single-kick evolution operator 

U = cxp {-ip 2 /2h) exp (-iKcos(q) fh) . (I) 

The expectation values are then computed from the time- 
evolved wavefunction. The value K c = 6.908745 • ■ ■ was 
found H to be a critical value in a wide parametric win- 
dow dominated by anomalous transport. We illustrate 
this regime by considering K* = 6.905 and contrast it 
with K = 11, where the effective diffusion coefficient is 
actually below the quasilinear value D q i (seen in Fig. |l|). 
The quantization scale h = 27ra, where several irrational 
values of a in the range 0.04 — 1 were considered. In 
the quasilinear limit, the quantum diffusion is expected 
to saturate at t* = aD q i/h 2 where a was numerically 
shown to be 1/2 ||l|,|lT||. £ measured in angular momen- 
tum quanta is ~ t*. 




FIG. 3. Variation of < p 2 > (normalized to the saturation 
value D q it*) with the number of kicks t. Two representative 
values of a, one large (Q,a ~ .61803399) and the other small 
(SC,o ~ 0.042340526), are shown for K = 11 (dashed lines) 
and K* — 6.905 (solid lines). Note that the larger a value was 
time-averaged to remove fluctuations that mask the trends. 
The inset shows a log-log plot (at longer times, w 5000 kicks) 
for a = 0.042340526, K* = 6.905 (upper curve) and K = 11.0 
(lower curve). The corresponding slopes are 0.21 (solid line) 
and 0.06 (dashed line) respectively. 

Figure. || shows < p 2 > /(D q it*) as a function of time 
for two representative values of a. For large a, both val- 
ues of K result in strong saturation with K = 11 having a 
higher saturation value than K* — 6.905. However, with 
decreasing a the enhanced diffusion for K* = 6.905 be- 
comes evident and there is a slowdown in growth in lieu 
of saturation (over the time considered). By contrast, 
K = 11 appears to saturate at a value which is about 
an order of magnitude lower. This difference is consider- 
ably larger than earlier expectations based on quasilin- 
ear analysis. The inset in Fig. displays the existence 
of power-law behavior in < p 2 > for longer times at K* 
as compared with saturation at K = 11. Note that the 
slope points to sub-diffusive rather than to super-diffusive 
growth for which we have no clear explanation at present. 
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A stronger signature appears on considering the vari- 
ation of localization length £ with K. It was shown Jll]] 
that £ tracks the variations in the classical diffusion 
coefficient D when plotted with respect to K q = 
Ksm(H/2)/(K/2). It is clear that this can be an im- 
portant correction for larger values of h. Figure ^ shows 
the variation in £, normalized to the quasilinear estimate, 
with respect to K q . £ is obtained by fitting an exponen- 
tial to the time-dependent wavefunction for long times. 
The peaks in the inset coincide with the classical accel- 
erator modes at K = 2n and 4-7T and £ computed at 
different times reflect the same features. The more de- 
tailed scan displays a reasonable correspondence with the 
results shown in Fig. [l], despite being far from the semi- 
classical limit. It should be noted that the highest value 
occurs at K q w K c defined earlier. 
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FIG. 4. Oscillations in the localization length £ nor- 
malized to the quasilinear estimate as a function of K q . 
a ~ 0.131267463 which makes the corrections to the classi- 
cal stochasticity parameter K significant. The inset shows 
peaks associated with the accelerator modes while the main 
figure zooms in on the region considered in Fig.l. The inset 
also shows £ computed after 3000 kicks (points) as well as 
9000 (line) kicks. Note that the wavefunction averaged over 
600 kicks is used to determine £. 

We now consider the nature of the localized wave- 
function when the classical transport is predominantly 
anomalous. Figure [5| shows the evolved momentum dis- 
tribution averaged over 600 kicks. Panel (a) corresponds 
to K q at a peak in Fig. ^ while (b) sits in the valley. Lo- 
calization does occur though the characteristic lineshape 
is strongly non-exponential for K q in the neighborhood 
of the peaks. The contrast with the clearly exponential 
lineshape shown in (b) is striking. It should be noted that 
the central region in (a) also displays curvature and an 
exponential fit results in a £ which is considerably larger 
than the quasilinear prediction. However, in case (b), 
the computed £ is well-approximated by the quasilinear 
estimate. Moving away from the peak in K q decreases 
the prominence of the shoulders in the lineshape though 
these are still clearly visible. We also note that the shoul- 
ders develop over a time (which depends on %) which is in 



the neighborhood of the cross-over time from exponential 
to power-law behavior shown in Fig. pi 
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FIG. 5. Lineshapes after 6000 kicks starting from 
a plane wave initial condition. (a) corresponds to 
K q m K c = 6.90845 ■ ■ ■ while (b) corresponds to K q = 10.69. 
a « 0.131267463 and the noise in the lineshapes was reduced 
by averaging over 600 kicks. 

Arguments for deviations of the localization length 
from the simple estimates proceed from the long- 
time averaged momentum distribution < f n >= 
\' l Pm(no)\ 2 \?Pm(n)\ 2 , starting from a plane wave n = 
uq initial condition. ip m (n) is the plane wave representa- 
tion of the eigenfunction with quasienergy oj m . The anal- 
ysis of localization assumes that all quasienergy states 
are exponentially localized and that the fluctuations in 
£ are small. If this were true in the case of anomalous 
transport, the oscillations in Fig. ^ would be directly cor- 
related with the number of quasienergy states participat- 
ing in the dynamics. To address this view, we computed 
that participation ratio PR from the time averaged au- 
tocorrelation function 

PR= lim i| < -0(T)[V(O) > | 2 . (2) 

1 — >oc ± 

As seen from Fig. |^, the only correlation with Fig. [| oc- 
curs in the valley where the lineshape is exponentially lo- 
calized. At the peaks, these results strongly support the 
idea that anomalous transport results in strong fluctua- 
tions in £ among the individual quasienergy states. We 
have verified this by constructing the 'entire' spectrum 
of the quantum map U, in the momentum representa- 
tion, under conditions necessary for dynamical localiza- 
tion. The average value of £ does not reflect the large 
changes in D though dispersion A£ exhibits large fluc- 
tuations. This view is being explored further at present. 
As we argued earlier, a similar scenario may apply to the 
interplay of scarring with dynamical localization. 
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FIG. 6. Participation ratio as a function of K q . Time av- 
erages of the autocorrelation function over lOOO(points) and 
6000(line) kicks are shown, a w 0.131267463. 

Theoretically, the classical flights of arbitrary length 
result (for the standard map) from a hierarchy of bound- 
ary layer islands |^. The island area S n , associated pe- 
riod T n and Lyapunov exponent r„ scale as: S n — XgSo, 
T n = A T Tb an d r„ = ApTo with the generation n of the 
island chain. As < 1, At > 1 and Ar ~ 1/Ay. It is 
straightforward to construct a semiclassical break-time 
appropriate to anomalous transport using these scaling 
rules. 

Quantum flights are possible as long as S n > % lead- 
ing to a critical value n% — ln?i/5o|/| In As|. The 
corresponding time is Tn — T Q (So/h) 1 ' ' M where fi = 
In As | /In At | > 1 is a classical transport exponent ||. 
We generalize the break time Tn defined earlier to 
ln(J/?i)/r n where / was some characteristic action. The 
longest time corresponds to the smallest r„ = Tn = 1/7}, 
from which we have 

r n = T h In (I/h) = T (So/h) 1 ^ In {I /ft) . (3) 

The result points to a power-law dependence on H rather 
than the logarithmic one as well as a possible crossover 
in scaling behavior. Note that the saturation time Tn 
for flights is a pure quantum effect and one that is nec- 
essary for localization, which cannot occur if arbitrarily 
long flights exist with high probability (not exponentially 
small). Theoretical estimates are considerably harder 
when fi > Sq. 

The issue of fluctuations in £ is relevant to an atom 
optics realization of the kicked pendulum Here, the 
localization properties are similar to those in the quan- 
tum standard map, with the complication that small is- 
lands persist in the classical dynamics. As we have il- 
lustrated, this could lead to strong fluctuations in the 
localized wavefunction. More recent experiments have 
realized the delta-kicked rotor system Jl3| where the mea- 
surement of diffusion, over timescales considerably longer 
than those predicted earlier, can provide direct evidence 
of anomalous transport in select windows of K. We note 
that the signatures of anomalous transport are similar 
to those resulting from external noise, in both the line- 
shape and the time evolution of the energy. This may 



be especially relevant in the weak noise limit and will be 
considered in detail elsewhere. 

In conclusion, we have presented evidence of the role of 
anomalous diffusion in enhancing quantum fluctuations. 
Note that the standard map is not ideal for isolating the 
effects of anomalous transport. The kicked harmonic 
oscillator, which results in the classical web map, is a 
more interesting candidate where super diffusion domi- 
nates the classical dynamics |(| . The absence of localiza- 
tion in that system Jl4| and a recent suggestion for an 
experimental realization in the context of ion traps |l~5| ] 
makes this an attractive direction to pursue. 

The work of BS was supported by the National Sci- 
ence Foundation and a grant from the City University of 
New York PSC-CUNY Research Award Program. GMZ 
was supported by U.S. Department of Navy Grants Nos. 
N00014-93-1-0218 and N00014-97-1-0426. The authors 
are grateful to Mark Edelman for assistance with the 
Poincare recurrence plot. 



[1] see, for example, Quantum, Chaos : between order and dis- 
order : a selection of papers, edited by G. Casati and B. V. 
Chirirov (Cambridge University Press, New York, 1995). 

[2] S. Fishman, D. R. Grempel, and R. E. Prange, Phys. Rev. 
Lett. 49, 509 (1982);D. R. Grempel, R. E. Prange, and S. 
Fishman, Phys. Rev. A 29, 1639 (1984). 

[3] G. P. Berman and G. M. Zaslavsky, Physica 91A, 450 
(1978). 

[4] E. J. Heller, Phys. Rev. Lett. 53, 1515 (1984). 

[5] E. J. Heller and S. Tomsovic, Physics Today, July 1993, 

p. 38 and references therein. 
[6] G. M. Zaslavsky, M. Edelman and B. A. Niyazov, Chaos 

7, 159 (1997). 

[7] S. Benkadda, S. Kassibrakis, R. B. White, and G. M. Za- 
slavsky, Phys. Rev. E 55, 4909 (1997). 
[8] R. Ishizaki, T. Horito, T. Kobayashi, and M. Mori, Progr. 

Theor. Phys. 85, 1013 (1991). 
[9] R. Roncaglia, L. Bonci, B. J. West, and P. Grigolini, Phys. 

Rev. E 51, 5524 (1995). 
[10] R. Ketzmerick, K. Kruse, S. Kraut, and T. Geisel, Phys. 

Rev. Lett. 79, 1959 (1997). 
[11] D. L. Shepelyansky, Physica D 28, 103 (1987). 
[12] J. C. Robinson, C. Bharucha, F. L. Moore, R. Jahnke, 

G. A. Georgakis, Q. Niu, M. G. Raizen and Bala Sun- 

daramRhys. Rev. Lett. 74, 2974 (1995). 
[13] F. L. Moore, J. C. Robinson, C. F. Bharucha, Bala Sun- 

daram, and M. G. Raizen, Phys. Rev. Lett. 75, 4598 

(1995). 

[14] G. P. Berman, V. Yu Rubaev, and G. M. Zaslavsky, Non- 
linearity 4, 543 (1991). 

[15] S. A. Gardiner, J. I. Cirac, and P. Zoller, Phys. Rev. Lett. 
79, 4790 (1997). 



4 



